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Motivated by their gravitational-wave driven instability, we investigate the influence of the crust 
on r-mode oscillations in a neutron star. Using a simplistic model of an elastic neutron star crust with 
constant shear modulus, we carry out an analytic calculation with the main objective of deriving 
an expression for the "slippage" between the core and the crust. Our analytic estimates support 
previous numerical results, and provide useful insights into the details of the problem. 



I. INTRODUCTION 

\b'- 

i The study of neutron star oscillation modes has long been a vivid research area, attracting additional attention 
' in recent years after the discovery that inertial modes of rotating fluid stars are prone to unstable behaviour when 
I coupled to gravitational radiation This instability is characterised by growth timescales short enough to make it 

astrophysically relevant. For example, unstable r-modes (the class of purely axial inertial modes) have been proposed 
' as a possible mechanism for limiting the spin-periods of neutron stars in Low-Mass X-ray Binaries (LMXBs) and 
as a potential source of gravitational radiation for the ground-based detectors. For a detailed review of the r-mode 
■"sj" ' instability we refer the reader to Ref. 

. Studies of the r-mode instability have revealed a host of possible damping mechanisms counteracting the growth of 
an unstable mode. Among the most efficient is the damping due to friction, via the formation of an Ekman boundary 
] layer at the crust-core interface, see for a recent discussion. For mature neutron stars, this mechanism could be 
the key damping agent for unstable r-modes. The crust forms once the star cools below 10^° K or so, soon after the 
■ star is born. At the temperatures relevant for a mature neutron star, the heavy constituents in the core (neutrons, 
t~ — . protons, hyperons) are expected to be in a superfluid/superconducting state. This leads to the suppression of, for 
T— I ' example, the hyperon bulk viscosity H, IE t;he emergence of mutual friction (see also Q for a recent discussion 
, of the shear viscosity in a superfluid neutron star core). In this paper we focus on the role of the crust. A detailed 
' discussion of the natural first approximation, where the crust is modelled as a rigid spherical "container" , and the 
<^ viscous Ekman layer that forms at the core-crust interface, can be found in 

• In the present paper we explore a different aspect of the problem, namely, we consider a more refined model with an 
Oh| elastic crust characterised by a finite shear modulus. A model along these lines was considered b y L ee & Strohmayer 
Q . (hereafter LS) in their numerical study of oscillating neutron stars (for some earlier studies see jlEl and |0|). Being 
?H ' elastic, the crust can sustain forced oscillations driven by (say) a core mode. In addition the crust has its own unique 
c/3 , dynamics, with a set of oscillation modes associated with the elasticity. The problem becomes very rich since one 
^ ' expects resonances/ avoided crossings to take place when a core-mode frequency approaches a crust-mode frequency 
|l2l |. for example as the star spins down. This picture is certainly much more complicated than one with the fluid 
simply "rubbing" against a rigid boundary wall. 

If we model the crust in a more sophisticated manner we should also expect the relevant r-mode damping timescales 
. to be altered. This was emphasised by Levin & Ushomirsky (hereafter LU) ^3 ^^'^ used the formalism of LS but 
considered the simpler case of a uniform density, slowly rotating spherical star, with a constant shear-modulus crust. 
They argue that the change in the r-mode damping timescale, compared to the solid crust case, can be quantified in 
terms of the relative motion between the crust and the core. Introducing the dimensionless "slippage" Sc — {v — vo)/vo, 
where v is the core fluid velocity immediately beneath the crust and vq is the corresponding velocity in a fluid star, 
LU argue that an approximate revised damping timescale is the standard Ekman result multiplied by the factor 
LU attacked the problem by direct numerical integration of the crust-core equations, and showed that the slippage 
is Sc ~ 0.05 — 0.1 in a typical case. Hence, they found that the r-mode damping timescale could be more than 100 
times slower than for a solid crust. For rotation rates such that the inertial r-mode frequencies are close to the various 
shear modes in the crust, the slippage factor can however become much larger (of order unity). These results were 
confirmed in a later study by Yoshida and Lee [l^. 

Our aim in this paper is to study the interaction between the oscillating fluid core and the crust using purely analyt- 
ical tools. Our decision to pursue the problem analytically was driven by a desire to obtain a better understanding of 
the intricate coupling to the crust, and the way that the flnal result scales with the different parameters. Although our 
chosen strategy leads to a significant amount of algebra, we believe that the end result justifies the effort. Certainly, 
the analytical relations we derive illustrate the underlying physics much clearer than a direct numerical calculation 
would. To achieve our aims, we make use of boundary-layer theory ^ combined with a truncated expansion in spher- 
ical harmonics for the perturbed fluid variables. In order to keep the calculation tractable, we consider a uniform 
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density, slowly rotating Newtonian star with a constant shear modulus crust. We study in detail how the crust is 
driven by a single r-mode in the core, and then compute the slippage and the induced modifications to the r-mode. 
Our formulation of the problem is very similar to that of LU. The main difference is that we assume that a purely 
axial r-mode in the core acts as the driving force, while LU included 0{fl^) (polar) corrections to their core mode {il 
is the rotation rate). Our study is also limited to a crust with a relatively small shear modulus (more like a jelly than 
a solid). In this sense our analysis is less general than that of LU 13], but we show that realistic estimates of the 
shear modulus support our assumptions. Where both calculations are valid, our results are in good agreement with 
those of LU (they are not identical due to the slightly different formulations of the problem). 

The paper is organised as follows. In Section |n] we briefly discuss the elastic character of the crust, providing some 
justification for our assumption of a constant shear modulus. The main scope of the paper is developed in Section IITTI 
where we discuss neutron stars with elastic crusts, and study the delicate coupling between the core-flow and the 
crust. The concluding Sections. IIVI fcl\n of the paper provide a discussion of the physical relevance of our results and 
a summary. In these sections we compare our results to previous work on this problem. We also point out interesting 
extensions of the present work. 

Aware of the somewhat heavy algebra involved in our analysis, we have made an effort to make the sections 
discussing the results self-contained. This should make the results accessible also for readers who are not prepared to 
bear with us as we work through the details of the core-crust problem. Many of the technical details of our calculation 
are located in the Appendices. 

II. MODELLING THE NEUTRON STAR CRUST 

We first discuss in some detail the properties of an elastic neutron star crust. This is the case that is relevant for 
most neutron stars since the crust is thought to form already at temperatures near 10^° K, possibly within minutes 
after the neutron star is born. Furthermore, the crust is not expected to be particularly rigid. Even though it supports 
shear stresses its nature is more like jelly, and one would expect it to take part in any global oscillation of the star. 
In addition, it will support distinct crust modes characterised by the speed of shear waves in the nuclear lattice. 
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FIG. 1: The shear modulus, based on Eq. Q and data for the equation of state of Douchin and Haensel fl^ . The vertical 
dotted line indicates the location of neutron drip, while the dashed lines indicate the region near the base of the crust where 
the nuclei undergo a series of changes in shape, cf. the discussion in [T^. 

In order to keep the problem analytically tractable, we will adopt a simple uniform density neutron star model. 
In addition, we will take the fluid to be inviscid since we are mainly interested in the mechanical, dissipation free, 
coupling between the r-mode and the crust. Furthermore, we assume that the crust is characterised by a uniform 
shear modulus fi. This assumption (which was also adopted by LU (13.] ) is somewhat unrealistic HI, but it greatly 
simplifies the analysis. A study of the non-isotropic case requires a numerical solution. On the other hand, realistic 
models for the crust suggest that the speed of shear waves is remarkably constant throughout the crust. We illustrate 
this in Figure n where we have combined the standard expression for the shear modulus of a bcc lattice 

^ = 0.3711^1^;^ (1) 
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where Z is the number of protons per nucleus and n^v is the number density of nuclei, with the data for the recent 
equation of state of Douchin and Haensel 0| . As is clear from this figure the ratio 

1,= ^^ IQie cmVs' (2) 
P 

only varies by a factor of a few in the crust. It should also be noted that [i does not vanish at the core-crust 
phase-transition. Given the results shown in Figure ^ it is clearly not too unreasonable to treat /i as a constant. 

The crust's rigidity can be quantified by means of the ratio (essentially the ratio between the speed of shear waves 
and a characteristic rotation velocity, or the relation between the elasticity and the Coriolis force), 

The crust behaves almost like a fluid when A ^ 1 and is essentially rigid in the opposite limit. In the present work 
we will only be concerned with the former limit. In this sense, our analysis is less general than that of LU who put 
no restrictions on the ratio Q. Introducing a characteristic rotation frequency 



^1 = 9^^ ^IR^ « 1020 cmVs^ (4) 



(such that the Kepler break up rotation rate is about O.Gfio) we see that we will have 

A=»10-(t)^ (^) 

This clearly shows that the A <C 1 limit is appropriate provided that we do not let ^ 0. In a sense, this illustrates 
the fact that our calculation is limited to rotating stars. Another reason for this is that we will assume that the core 
fluid undergoes an inertial r-mode oscillation. This would not make sense for a non-rotating star as these modes then 
become trivial. 



III. CRUST-CORE COUPLING 



A. Formulating the problem 

The physical situation we wish to examine is that of a single core r-mode coupled to the crust. As a result 
of this coupling the crust will oscillate. At the same time there will be some "back-reaction" on the core mode. 
Mathematically, one needs to solve the Euler equation in the core region, and the corresponding equations for an 
elastic body in the crust region. Boundary conditions need to be imposed at the core-crust interface and at the 
star's surface and centre. Our aim is to approach the problem through a two-parameter expansion: a slow rotation 
expansion in together with an expansion in A. The discussion in the previous Section illustrates the reason why 

this is a tricky problem. Yet we will see that we can make signiflcant (possibly surprising) progress by combining 
these expansions with boundary layer techniques. 

We assume a uniform density, incompressible fluid. Since we work at leading order in rotation we can assume 
our star to be spherical. Moreover, we adopt the so-called Cowling approximation where the perturbation in the 
gravitational potential is ignored. This is a sensible assumption, since for the quasi-axial modes we are interested 
in the associated density perturbation (and consequently the perturbation in the gravitational potential) is a higher 
order correction. In the rotating frame, the linearised Euler and continuity equations for the crust flow are 

dU + 2nxdt£,^~-V5p + p,V^S, (6) 
P 

V • e = (7) 

where we have taken jl = fj,/ p = constant. For the present problem it is natural to work with the displacement vector 
5. The main reason for this is that the shear introduces a Hookean restoring force into the Euler equations, i.e. terms 
~ /iV^^ p. [Tfil[lll |. The core flow is described by the /i ^ limit of the above equations. 
We assume a decomposition of the form. 



r 

I > m 



sin 6 J \ sin I 



(8) 
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where the relative phase between the polar {Wi, Vi} and axial components {Ui} is chosen such that the final pertur- 
bation equations are real valued. We also decompose the perturbed pressure as 



^ l>m 



(9) 



The three components of the Euler equation ((BJ then reduce to (in the following, a prime denotes a radial derivative); 



{rh'i ~ a^Wi + 2manVi] Y," + 2(jWi smOdeYl^ + /I <^ 2^ - W^' + + l 



Wi 



Wi-2Vi 



= 



(10) 



5^ = ^ {-ma^ViY{^ + {hi - a'^Vi) ^^iiBdeYr" + 2mcrf7Vi coseF;™ + ^aWi cQsQfA^BdeYi 



mf, |[//' - 1(1 + 1)^1 Yr - A - 1(1 + 1)^ + 2 



— >>sm09ei^ 



(11) 



S"^ = [{"^^' + 2crrJWl - ma'^Vi} i^™ - ct^J/; sinOdeY/" + 2maWi cosflY;™ + 2(jnVi cosOsinedgYf 



v. 



2anWi cos^ OYl^ - mfl \ ¥{' - l{l + 1)-^ + 2-^| Y;'" - /i <{ [//' - l[l + 1)-| [> sin6'9eY, 



(12) 



These equations are identical to thosegiven in [olll ij. specialised for constant shear modulus and incompressible flow. 

As discussed by Lee & Strohmayer [9| it is convenient to consider two particular combinations of the above equations. 
The first is the radial component of the vorticity, 



and the second is the "divergence equation" , 

From cqn. (|13|l and from the combination of cqns. p4|l and (|10|l we arrive at the following equations, 



(7 

2£L 



/(/ + l)--,n 



C/i + ^ [ (/ + \)Qi{Wi-Y -{I- m-i} - lQi+i{Wi+i + {1 + 2)Vi+i}] 



= 



(13) 
(14) 



(15) 



— [/(/ + 1){{1 - l)QiUi-i -{1 + 2)Qi+iUi+i} -{1- 1)(/ + l)QirUl_^ - l{l + 2)Qi+irU'i+^ - m{rV{ - Wi}] 



rvr 



2Wl' + {l-m + 2)^ 



= 



where 



' {l + m){l - ra) 
{21- 1)(2Z + 1) 



(16) 



(17) 



The continuity equation jT)) provides one additional expression 

(rWiY - l{l + l)Vi ^ (18) 

Equations p5|l . (|16|1 and (|18|) are the basic equations for our problem. They involve only the eigenfunctions of ^ 
and the frequency a and are equivalent to the initial Euler equations. It is interesting to note that our equations 
only couple multipoles of order I to those corresponding to Z ± 1. This is in contrast to the (fluid) equations used by 
eg. Lockitch and Friedman 15] to study inertial modes, where coupling to Z ± 2 multipoles was also present. In this 
sense, the present set of equations is simpler and could possibly be advantageous also in studies of the fluid inertial 
mode problem. 
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B. Boundary conditions 

The dynamical equations are accompanied by boundary conditions at the crust-core interface (at r = Rc) and at 
the surface of the star {r — R). These boundary conditions can be formulated in terms of the traction components 

Tr = 2ndrC-^P (19) 



Tg = ^ 



r r 



r r sm t) 



(20) 



(21) 



where Ap is the Lagrangian variation in the pressure. The traction vector is constructed by projecting the fluid 
stress tensor along the unit normal vector at the interface, and represents the force per unit area. Perturbing (in the 
Lagrangian sense) this projection results in the expressions above. 

At the crust-core interface all traction components as well as the radial displacement vector need to be continuous. 
Specifically, this requirement leads to (after setting I — m in the axial component) 

(i) Vanishing tangential traction at the interface and surface: 

/i[rC/;, - 2U,n] =0 at r = R^, R (22) 

A ['^C+i + - 2Kn+i] =0 at r = i?„ i? (23) 

(ii) Continuity of the radial traction at the interface: 

5p''°'° = - 2ndrC ^ - - 2fL {^^^ at r = Rc (24) 

(iii) Vanishing of the radial traction at the surface: 

6p + rP ~ 2ndre = => R^hrn+i^-goRWm+i+2^i[RW;^+^~Wm+i] at r = R (25) 

where we have defined go = GM/R^, the surface gravitational acceleration (which appears after using the background 
Euler equation Vp = — pV$). One can show that this final condition reduces to the standard condition of vanishing 
Lagrangian pressure perturbation at the surface in the limit /i — > 0. 

(iv) Continuity of the radial displacement at the interface: 

Wrn+liR;)^Wm+liRt) (26) 

In addition, we require all core flow variables to be finite at r = 0. 

C. Simple mode solutions? 

Before we begin our main analysis, it is worth discussing the nature of some comparatively simple solutions. First 
of all, we ought to be able to reproduce the standard fluid r-mode solution from equations ((T5l) and l(TB|) . Let jl = 
and assume that the only non- vanishing displacement is U„i. Then eqn. (|15|l with I = m leads to 

[m(m + 1)0-2 _ 2m<Tn] [/„, = (27) 

and we must clearly have either = or 

2n 



m + 1 



(28) 



which is the standard result. Meanwhile eqn. (|16|l with / = to — 1 is trivially satisfied since Qm — 0. Finally we need 
eqn lfTC|l with / = to -|- 1; 

rU^-{m + l)Um = ^ Um = Ar"'+^ (29) 
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which is the expected eigenf unction. 

Given the perturbation equations from the previous section it is also easy to prove that inertial modes in uniform 
density stars cannot have polar terms as their highest multipole contribution Consider a mode which truncates 
after I — n (say) and for which the n**^ terms are polar. Then eqn. 115|l with I = n + 1 provides the relation 

Wn -nVn = (30) 

Combining this with the continuity equation we see that we must have 

= r" (31) 

which cannot satisfy the required surface boundary condition Wn{R) = 0. In other words, the highest multipole term 
of an inertial-mode solution cannot be polar. 

More important for the present investigation is to ask whether we can have purely axial modes in a star with a 
fluid core and an elastic crust. We have already shown that in order to have only Um (say) nonzero in the fluid core, 
we must have a mode with the r-mode frequency. We also have the eigenfunction from Eq. (|15|l . In order to be 
consistent, this eigenfunction must satisfy the remaining part of Eq. (|15|l in the crust, i.e. 

f/;;;-m(TO + i)-^ = (32) 

This is, in fact, the case. But our solution is nevertheless not acceptable because it cannot satisfy the boundary 
conditions of vanishing traction at the base of the crust and the surface of the star [Eq. (|22|) above] . We require 

Ul^ - 2— = at r = i?c and r = R (33) 
r 

which clearly is incompatible with our solution. Thus we conclude that a purely axial mode is no longer possible in 
the presence of the crust. 

This is, in fact, a bit puzzling because we can extend the argument to demonstrate that a crust mode with only 
[C/jp, Kn+i, Wm+i] 7^ is also inconsistent! This would seem to cast some doubt on the numerical results of LU 
|l3| who onlyinclude these terms. In reality the situation is a bit more complex. As also mentioned by Lee and 
Strohmayer [3, the problem should be approached in a way similar to the inertial-mode problem for a fluid star 
[T5I I. Even though one should not at the outset restrict the number of niultipoles included in the calculation one 
will always in practise truncate after some multipole. Since this involves neglecting some of the relations implied by 
the recurrence relations it renders the solution somewhat inconsistent (as in the case of LU). However, one would 
expect the mode-solution to converge when more multipoles are accounted for, and it does not seem unreasonable to 
expect that a lower order representation will be fairly accurate. Indeed, this turns out to be the case in a similar 
calculation for the standard Ekman layer problem where the known exact solution is quite well approximated by a 
[Um, Vm+1, Wm+i] truncated representation (see Q for discussion). As we are seeking an analytical approximation to 
the solution, it is natural to truncate the solution at the simplest "meaningful" level. This is an obvious weakness of 
our analysis. A more detailed study, investigating for example the convergence as more multipoles are accounted for 
and whether the outcome depends on the parameter values, requires a numerical solution. 



D. Core flow 



As we have mentioned, we will assume that the bulk core flow is described by a single I — m r-mode, which is 
capable of penetrating the crust region and driving the flow there. In return, it picks up corrections due to the coupling 
with the crust. Furthermore, we will truncate the multipolar expansion (see can. ( 15 ) )) at the level {Um, Vm+i, Wm+i}- 
Note that the same truncation was adopted in previous studies of the problem [3, ■ 

For the dynamical variables we assume expansions of the form 

Q = Q° + Q' + Q^ + ... (34) 

where, for brevity, we have collectively denoted Q = {Um.,Vm+i,Wm+i,hm+i,cr}. Expansion H34() has no obvious 
small parameter in it; it is written anticipating the fact that there is a one to one correspondence with the similar 
expansion in the crust region (which has the well defined expansion parameter /i, see Section Fill Ell . The reason for 
including terms up to second order will become apparent in the discussion of the crust flow. 

As we are assuming an r-mode solution at leading order we take the leading order polar pieces to be zero, i.e. 



(35) 
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The non-zero axial piece is normalised in such a way that A = 1 in (23), i-e. = at the base of the crust. 

The next step is to insert the expansion (|34|l in the main equations (|15|) - (|18|l . The resulting first and second order 
equations are listed in Appendix The analysis in that Appendix provides the results, 



m+l 



(36) 



(m + l)T/„\+i 



(37) 



m+l 



and 



[(72 + l<Jl] 

2n 



(to + l)Qm+ir 



m+l 



W 



Ui = Sr 



m+l 



(38) 

(39) 
(40) 
(41) 



where 7 and S are constants. 

These results demonstrate the importance of the frequency corrections cti and (T2 . Had we neglected them, we would 
have effectively killed all corrections to the core flow. Hence, the secondary core flow is linked to the crust-induced 
correction to the mode oscillation frequency. 

Finally, we need to calculate the pressure correction, anticipating the fact that the pressure appears explicitly in 
some of the boundary conditions. For this purpose it is convenient to use the £^ component of the Euler equation, 
eqn. ((TT|l . We easily obtain. 



''■m+l 



"m+l 



''m+l 



_ mQm + 1 rn+1 



(to + 1)2 
?m+l 



(to + 1)2 
m+l 



|^(to + 1)(2to-1) 



- TO7 



„m + l 



(to + If 



(to + 1)(2to, - 1) 



0'2 + 7^1 



+ m5 + (to - 2)(to + 1)^ ( ) 



„m+l 



(42) 
(43) 
(44) 



In the end, all core-flow eigenfunctions (up to second order) are expressed in terms of the (as yet unknown) parameters 
fij f2, 7 and 5. As we will see, these are linked to the crust flow through the interface conditions. 



E. Crust flow: setting up the equations 



In the crust region the Euler equations have additional term of form jiS/^S^. These terms are similar to the viscous 
terms in the Ekman layer problem but here they represent the crust's elastic restoring force. Despite the additional 
terms being dissipative in one case and oscillatory in the other, the problems are mathematically similar. The crust 
equations are likely to exhibit "boundary layer behaviour" when the coefficient /i, which multiplies the highest order 
spatial derivatives, is treated as a small parameter (i.e. in the limit A ^ 1). This feature will turn out to be of key 
importance when solving the equations. 

Since we are only considering an "almost fluid" crust, one would expect the crust flow to be primarily axial in nature, 
as it is driven by the quasi-axial core oscillation. In the crust, as in the core, we assume from the outset that Vm+i and 
Wm+i vanish at leading order, and truncate the multipole expansion © by including only the {Um,Vm+i,Wm+i} 
eigenfunctions. 

In the crust, the continuity equation H18|) with I = m + 1, gives 



{rW^+iY = (to + 1)(to + 2)Kx+i 
Meanwhile H15() with / = to, and (|16|l with I = m + 1, can be written 



(45) 



2n 



Urr 



U':,-m{m+l)^ 



— g,„+i + (to + 2)y„+i] = (46) 



8 



[(m + l)(m + 2)Q™+i[/™ - (m + 2)g,„+irC/;, - rV^,+i + W, 



m+l I 



rt/„+i - 2VK„+i + m(m + 3) 



m+l 



We now adopt a standard uniform expansion 17], 

g = g° + Q" + /i"{g' + g'} + + Q'} + 



(47) 



(48) 



where we have again used the collective notation Q = {UrmVm+i,Wm+i,hm+i,<y}- "Boundary layer" variables are 
denoted by tildes (and for them we assume dominant radial derivatives dr ~ 0{jl~°')). The remaining quantities 
pertain to the smoothly varying, background crust flow. Note that in 14811 the power-law exponent a is initially 
undetermined. Its value is fixed by the equations themselves, in order for them to make sense. The analysis of 
Appendix IbI shows that a — 1/2. 

Inserting the expansion (|48|) in eqns. H45|l - (|47|l we find two families of equations, one for the "boundary layer" and 
one for background variables. We derive all the dynamical equations required for the determination of the crust flow 
up to 0{jl) (the reason for pursuing the calculation to this order, rather than simply to ©(/i^/^) will become apparent 
below). Next, we expand the various boundary conditions (|22ll - (|2t)() . In order to avoid overloading the reader with 
algebra we carry out these calculations in Appendices El and 



1. Leading-order flow 

Having written down all the relevant equations, we can now solve eqns. (IBlOp and (|B15|) for the leading order flow. 
After decoupling these equations we end up with, 

~2y^n„ _ 160Sn(m + 2) , _ 

'"+i+(m + l)2(m + 2)'^™+i (m+l)4 ^^^^ 

Trying the Ansatz, V^^i ~ e^^-^i) we find for the exponent a pair of real and a pair of imaginary roots, 

x„.±^( (50, 

/i^/^ \ (m + 2)(m + 1)^ / 
''Jl^ Um + 2)(m+ 1)^ 

a=l + m{m + 2fQl^^^ (52) 



«A3,4 = ±^ . , OA/ , (51) 



where 



The occurrence of purely real and imaginary roots deserves further discussion. The former will lead to a pair 
of growing/decaying exponentials while the latter implies a pair of oscillatory exponentials. This is not the typical 
situation in boundary layer problems. Recall, for example, that in the Ekman problem Q we have only decaying and 
growing exponentials. The difference in the present problem is associated with the fact that the /iV^^ term in the 
Euler equation represents an elastic restoring force, not a viscous damping force. Intuitively, the purely oscillatory 
solutions to represent the elastic degree of freedom of the crust, eg. allows for the presence of the crustal shear 
modes in the calculation. This means that in order to properly construct the complete flow we need to keep all four 
solutions. Therefore we write, 

y^+i = A^e^^^ + B^e-^^^ + cos(A3C) - sin(A3C) (53) 

where ^ ^ r ^ Re- Similarly, 

U^^ = Auc^'^ + B„e-^i? + C„ cos(A3C) " MhC) (54) 



9 



As a result of eqns. (jB10|) and (|B15|) . the coefRcients are coupled in the following manner, 

fiXlim + If I = An^m + 2)Qm+i | | (55) 



flXiim + ly I I = 4f^^(m + 2)Q„,+i | I (56) 

We can now impose the boundary conditions (|Cip and ljC4|l at r = i?c- These translate into 

Xi{A„-B„) = (57) 

Ai(A„-B„) = A3-D„ (58) 

It is easy to see that these two relations are incompatible with (|55|l and H5t)|) . The only remaining option is the trivial 
result Du = Dy = 0, which in turn means that 

Au = Bu, and A„ = (59) 

Enforcing the conditions IjCip and ljC4|) again, this time at r = i?, leads to 

2AiA„sinh(AiA) = A3C„sin(A3A) (60) 

2AiA„sinh(AiA) ^ A3C„sin(A3A) (61) 

where A = R — Rc is the crust thickness. Combining these relations with the previous ones we are led to the trivial 
result, 

C, - C„ - A„ ^Ay^Bu = By=0 (62) 

However, the argument that non-trivial solutions do not exist is only true provided that 

sin(A3A) 7^ => A3A 7^ nn, n integer (63) 

This condition is obviously violated when il ^ 0. However, as we have already emphasised, we are not allowed to 

approach the nonrotating limit in the present framework (since for a given fi we always demand that A ^ 1). On 
the other hand, non-zero values of 51 that violate H63|l could well be accessible. Later on in this paper we discuss the 
physical meaning of these particular frequencies. 
We conclude that when (|63|l holds we must have 

U^m = = (64) 



It also follows that h^^i = 0. In fact, by using eqns. (|B10|) and ljB15p in (jB20|) . we can prove that h^j^i — 
identically. Hence, somewhat unexpectedly, we conclude that at the anticipated leading order all eigenfunctions must 
vanish. 

Before we move on, we should perhaps emphasise why we think it was worth working through the above calculation 
even though the end result was trivial. Basically, we wanted to highlight the exception H63(l . i.e., the existence of 
non-trivial solutions for certain parameter values. It is also worth pointing out that we would not necessarily have 
expected that the solution would be trivial at this level of approximation. Having said this, the result could have been 
inferred from the fact that the corresponding boundary conditions involve only "boundary layer" quantities [l7j . 



2. First order crust flow 

Having obtained a trivial result for the leading order flow, we must proceed to the next level (i.e. Oi^fj}/"^)). Let us 
first focus on the "boundary layer" variables. The continuity equation (|B3p gives, 

WUi = (65) 

We can also calculate V^j^i and U}^ from the radial vorticity equation IjBlip and the divergence equation (|B16|) . In 
view of the previous results, this last equation takes the form of eqns. (|B10|) and l|B15|l with the simple changes 
f/^ IJ^ and V^i^i — > The coefficients are (again) related by eqns. (|^ and ifS^ . 
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The relevant boundary conditions are equations ljC2|l and (jC5|l at r 

(to - l)AA2i?™ 



B„ 



A, 



2/21/2 sinh(AiA)Ai(A2 + A2) 

(to - l)A\lW^ 
" 2/21/2 sinh(AiA)Ai(A2 + A2) 

(to - l)AA2i?'" 
/ii/2sin(A3A)A3(A2+A2) 

(TO-l)AAf 



i?c and r = R. Their combination leads to, 



§) cos(A3A) 



/2i/2A3(A? + Ai)' 

Note that in the limit of a thick crust we have e'^^^ ^ 1, in which case we get Am Ay 



(66) 
(67) 
(68) 
(69) 



i.e. the growing 

exponentials are eliminated. Analogously, the two solutions that grow away from the surface will also vanish in this 
limit. 

The remaining components V^^j^i^ U}^-^ and W^j^i can be obtained from equations (|ij6|l . (|ijl2() and (|ijl7l) . Neglecting 
all homogeneous solutions, since we are only interested in the flow driven by the leading order r-mode, we find 



0n = (70) 
^„\+i - ^^(m+l)Q,„+ir™+i (71) 

W^^+i = {m + l)V,i+, (72) 



For the pressure eigenfunctions we again find that = identically (due to eqns. IjBlip . ljB16p and the vanishing 

leading order flow). On the other hand, 

1 2»AQ^+ig-i (2to2 + TO - 1) 

^n+l- -1/2 (^+1)2 ^ ^'-^l 

The flows in the neutron star core and crust "communicate" via the interface conditions. In view of our previous 
results, some of these conditions reduce to trivial identities. The remaining non-trivial conditions are eqns. (jC12|) and 
"1311 which flx A and 7, 



A = l, 7 = (74) 

The requirement for the radial traction to vanish on the surface gives, at leading order, 

C+i(i?)=0 (75) 

This familiar condition cannot be satisfled by our leading order solution for the simple reason that we have assumed a 
purely axial r-mode background flow in the core. Indeed, in order to be fully consistent, one would need to introduce 
0(^2) polar corrections, as for example in These pieces are neglected in our calculation and consequently the 
above surface condition should not be taken into account. 
At the next order we have. 



.90 

"■m+l — 

from which we obtain an expression for ci 



hl+, = -^W^^+i (76) 



An^ai{2m^ +m~l) = —(7i{m + lf (77) 
R 

The only way to satisfy this condition is by making cti = 0. Hence, the leading order correction to the frequency 
appears at 0(/2) (or higher). This result has the consequence that for the induced core flow, 

Kn+i - ^U^^O (78) 

At this point we have used all available dynamical equations and boundary conditions at ©(/t^/^). Our results 
allow us to compute the crust-core slippage, a crucial quantity for the determination of the Ekman viscous dissipation 
rate associated with the presence of the crust |l3j . Unfortunately, as we have just found, the crust-induced frequency 
correction is undetermined at this level. If this quantity is desired, one has to extend the calculation to the next order. 



11 



3. Second-order crust flow 



We have already worked out the second order corrections to the core flow. In this section we deal with the 
corresponding crust corrections. To begin with, from the continuity equation ljB7|) we immediately obtain, 



(79) 



where we have defined the integrals, 



h{r) 
hir) 



dr' 



Ra ^ 

dr 



hir) 



dr^ 



- cos[A3(r' - R,)] h{r) = / — sin[A3(r' - R,)] 



,-Ai(r'-_R,) 



dr' 



(80) 
(81) 



and we have used the continuity of the radial displacement ljC9|l to fix the integration constant such that Wm+ii^c) = 
0. Given this, condition ljC12() requires that 



/i^+i(i?,) = 

From eqns. (|B8|) . ljB13|l and ljB19|l we get (again neglecting all homogeneous solutions). 



2 _ (™+ l)Qm+l^ ^+1 



w 



m+l 



(m+ 1)^,2+1 



(82) 

(83) 

(84) 
(85) 



Moreover, from eqns. (|B23|) . ljB24|l in combination with eqns. IjBlSII . ljB7|) and ljB14p we obtain the pressure 
corrections, 



ft2 



-1 , ,x~i/o [{Qm+i^« + AJ/i(r) + {Q.,n+iBu + B^}l2{r) + {Q™+iC„ + C^}h{r) 

(m + Ijfj,'-'^ 

{Q,n+iDu + D,}h{r)] (86) 



2ilQ„i+i 2 - m(2m + 3) 
/i (m + l)(m + 2)'^^'' 



(87) 



The only remaining non-trivial condition (at this level of approximation) is equation (|C15|I which can be solved to 
provide the frequency correction 



0-2 



~^-(m+l) 

2nQ 



m+l 



2 - m(2TO + 3) 1 



(m + l)(m + 2) 4 



-I -1 



IV. FINAL RESULTS 



After a somewhat involved calculation, we are now in a position where we can summarise the results and extract 
physically interesting quantities. The two most relevant quantities are i) the core-crust slippage, which enters in the 
estimates of the Ekman layer damping rate of core r-mode oscillations, and ii) the crust-induced frequency correction. 
The latter is interesting because it will unveil the presence of parameter values where our approximation scheme 
breaks down. As we will discuss, this happens because of the existence of distinct shear modes in the crust. 

We begin by considering the extent to which the crust moves along with the core fluid. From the studies of the 
Ekman layer problem (see jj for references) we know that the potential presence of a discontinuity at the core-crust 
interface could have significant impact on the viscous damping of a core oscillation. In principle, the role of viscosity 
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r/R 

FIG. 2: Eigenfunctions Um and Vm+i as functions ol r/R for the Z = m = 2 r-mode and for Q, = CSJlo, jl = 10~*-R^no. 
Dashed and soUd curves represent the background and 0{ji'^'^) corrected eigenfunctions, respectively. The vertical dotted line 
represents the crust-core boundary, here taken to be at Rc = 0.9/?. 




is to smooth out any discontinuities in the inviscid solution. To date, there have been no detailed studies of the 
viscous problem for a star with an elastic crust. Instead the effect that the elasticity of the crust has on the Ekman 
damping timcscale has been estimated by the core-crust "slippage" ^3 • The magnitude of the relative displacement 
between the core and the crust is (obviously) given by 



(89) 



The relative contribution of the polar and axial pieces can be seen in Fig.Elwhere we show the m = 2 eigenfunctions 
for a fixed angular frequency = O.3f2o (about half the Kepler break up rate) and jl = 10~^i?^r2Q. For this particular 
rotational frequency the axial piece is dominant, although we note that the polar piece becomes dominant once we 
increase the rotation rate and cross the first of the crust-core resonances discussed later. It is evident from the results 
shown in Fig. El that the corrections to the basic I = m r-mode oscillation are everywhere small, thus providing 
credence to our approximation of truncating the multipole expansion. 

The figure also clearly illustrates the discontinuity at r = Re- We define the crust-core slippage Sc by normalising 
the corresponding relative displacement with respect to the standard r-mode displacement ^o- This leads to 



Co 



1 



m+l 



I Um {Rt )-U^iR-)\ +1 Vm+1 (Rt ) - Vrn+l {R' ) I 



1/2 



(90) 



For a rigid crust the slippage is obviously unity, since A^ = if the fluid satisfies a no-slip condition at the base 
of the crust. Within the framework of the present calculation (a crust with A <C 1), we can use the results for the 
0{fi}-/'^) flow fSections IIII El fc IIII D|) and write the leading order expression for the slippage. We have found that the 
corrected global r-mode eigenfunctions are, 



Um{r) 



6(0 p}'^ [A„e^^« -f B„e-^^« + C„cos(A3C) - D^sin{\^C)\ 



(91) 



Vm+i{r) = 6(0 fi^'^ [A,e^^^ + B., 



cos(A3C) - sin(A3C)] 



(92) 



Wm+i{r) = + O(/i) 



(93) 



where Ai and A3 are defined by (|50|) and H51|) . and the coefficients Au etcetera follow from H66|l - (|69() and H55|) . H5t)|l . 
6(0 represents the unit step function and as before Q — r — Re- From the above formulae we see that the discontinuity 
at the core-crust interface follows from 



UmiRt)-Um{R-) - il'/^ir^iRc) = fi^/^A^ + Bu + Cu) 

Vm+i{Rt)-Vm+i{R-) = Mi/'v;\+i(i?e) = M'/'(v4„ + B„ + a) 



(94) 
(95) 
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Then 1)90(1 becomes, 



5,= 



m+1 



1/2 



(96) 



The way the shppage enters the r-mode problem has been discussed by LU |0| and Kinney & Mendell . They find 
that the shppage enters quadraticaUy in the r-mode damping formula, which leads to a revised Ekman layer timescale 



TE 



Te 



(97) 



In Fig. 121 we show the slippage as a function of il, for m — 2 and for two values of the crust thickness. The results 
show that a typical value is Sc ~ 0.05. In other words, we estimate that the r-mode damping time is about a factor 
400 longer in the case of an elastic crust that partakes in the mode oscillation than in the case of a solid crust. This 
estimate is in good agreement with that obtained by LU. 




0.3 0.4 



FIG. 3: The slippage Sc as function of Cl/Qo for Rc — 0.9R (top) and Rc = 0.8R (bottom). Note that the shear modulus p, is 
expressed in units of R'^Qq. The vertical dotted line marks the Kepler break-up limit Qk ~ O.Gfio- In both panels we have set 
I = m = 2 for the core r-mode oscillation. 

Fig. O also shows the presence of resonances at certain discrete values of the angular frequency. These peaks 
correspond to, cf. (|^ . 



A3(r2„)A — TiTT, n integer 
Solving this condition we find the resonant frequencies, 



„ -1/2 '™ / , 1 N + ^ 



(98) 



(99) 



14 



These resonance points are apparent if we graph the frequency correction from eqn. H88|) . see Fig.0J 

Clearly, the number of resonant points in the interval f2 £ [0, fio] increases with increasing crust thickness and 
decreasing /i. In the vicinity of the resonances our calculation breaks down. This means that, in the vicinity of these 
parameter values our starting assumption that the motion is dominated by the Coriolis force and can be represented 
as an r-mode plus a small correction is no longer valid. One would expect this to happen when the core r-mode 
frequency is close to the various crust mode frequencies. In other words, the resonances in Figs. |2| and 01 indicate the 
occurrence of avoided crossings between core and crust modes. That this interpretation is correct is clear from the 
results of LU ,1^] and Lee & Yoshida J^] . Away from these points the mode frequency is essentially that of the core 
r-mode, in agreement with the numerical results of LU. Having said this, the careful reader may note that our small 
frequency correction differs in sign from that of LU for fast rotation. That the two results differ in this regime is not 
too surprising, and is likely due to LU including some higher order rotational corrections. 
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0.02 
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-0.04 
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FIG. 4: Mode frequency correction 02 (normalised in units of the leading order frequency cro) as function of r2/f2o, for I = m = 2 
and Rc = 0.8R. The displayed value for jl is in units of R^Qq and the vertical dotted line marks the Kepler break-up limit. 



V. CONCLUDING REMARKS 



In this paper we have presented an analytic investigation into the interaction between an r-mode core oscillation and 
an idealised elastic crust for a uniform density, incompressible, rotating neutron star. The essence of our calculation 
is the use of boundary layer theory techniques. Our results are in agreement with previous numerical treatments 
of the problem [H El. They also offer a more detailed description of the workings of the crust-core coupling and 
the emergence of resonances between the r-mode and the various crustal modes, which leads to avoided crossings in 
the numerical studies (l2l Il3j . We have provided simple analytic results for quantities like the crust-core slippage 
and resonant frequencies (eqns. (|9ti|l and (|99|l ) which are of importance if one wants to estimate the viscous damping 
timescale for an unstable r-mode. 

There are several ways that one can refine and/or expand the present analysis. Recall that we have taken the fluid 
to be inviscid. This assumption could be relaxed and the resulting formulation should allow a consistent computation 
of the viscous dissipation rate. Note that presently this dissipation has only been estimated by rescaling the Ekman 
damping rate with the slippage factor |l3| . Other ways to improve our model would be to consider a crust with a 
non-uniform shear modulus (along the lines of '^), and/or a non-spherical rotating star. Both options would make 
an already complicated calculation even more complicated, and one would have to approach the problem numerically. 

An interesting extension would be the description of a neutron star with an elastic crust and a magnetic field. Based 
on the results from the study of the magnetic Ekman problem T^jT^ it is clear that a magnetic field may significantly 
change the associated r-mode damping timescale. The combination of an elastic crust and a magnetic field could well 
provide the most important damping mechanism for any unstable stellar modes. So far this combination has been 
investigated in 'l8l| , but not within a fully consistent framework. The crustal elasticity was only taken into account in 
the form of a slippage factor rescaling (as in (|97|l 'l the magneto- viscous dissipation rate computed for a solid boundary. 
Moreover, the core oscillation is assumed to be that of a normal, non-magnetic r-mode, even in the presence of strong 
magnetic fields which would most likely change the mode structure. There are also difficult issues of principle since 
the interior magnetic field structure is essentially unknown for neutron stars. 
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One would also have to account for the potential presence of superconducting protons in the core. A closely related 
problem arises from the fact that the crust is expected to be penetrated by superfluid neutrons. So far there have 
been no attempts to account for this feature, even though one would suspect that it could have both quantitative and 
qualitative impact. 

From this list of outstanding problems it is clear that, after quite a lot of algebra, we have only made some small 
progress towards the understanding of an area where many serious challenges remain. 

APPENDIX A: CORE FLOW 

Using the expansion in the basic equations ((T5|l . and l(TH|l. when applied in the core (i.e = 0), we obtain 
at first-order, 

(m + l)ai[/0 = 2f]g„,+i[W-,i,+i + (m + 2)l-„\+i] (Al) 

m(m + 2)Q„,+i [rU^ - {m + 1)U^^] = 2rC+i - mrW^+, - {m + 2)W^^, + m{m + l)(m + 2)V^^, (A2) 

{rWl+iY = (™+l)(m + 2)T4+i (A3) 
The corresponding second-order equations are, 

(m + 1)^72(7° +(m + l)ai[/^ = 2nQ„^+i[Wl+, + im + 2)V^^,] (A4) 

m{m + 2)Q„,+^[rU^ - (m + l)C/2„] - ^(rn + l)(m + 2)[rV,l:+, ~ W^^^] + m(m + l)(m + 2)Vl^^ 

-h 2rV;2Vi - mrW^^Vi " ("^ + 2)M^™+i (A5) 

From eqn. IjAip we get, after using the continuity equation H18() and disregarding the homogeneous solution which 
becomes singular at the origin, 

W^™+i = |^("^ + lfQ™+i^-™+' (A6) 



Then eqn. (fT^ gives. 



Ki+i = |^("^ + l)Q™+ir"+i (A7) 



Solving the remaining first-order equation l|A2p we get, 

Ul = jr^+' (A8) 



where 7 is a constant (since the source term in (jA2|l vanishes identically this is the homogeneous solution) . 
Similarly we find for the second-order flow, 

Wl^, = ^^^^i^g^[..+7-i]r™+^ (A9) 
= ^-^^^^^^^W2+ia,]r-+^ (AlO) 
Ul = ,5r'»+i (All) 

with S another integration constant. 

APPENDIX B: CRUST FLOW 

In this Appendix we insert the expansion (|48|) in the main dynamical equations (|46|l . H47|) and H45|) . We group 
together terms of the same order in fl and then split those groups into "boundary layer" and background pieces. For 
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clarity, we will discuss this process in detail only for the continuity equation (which is the simplest of our equations). 
For the remaining dynamical equations we just list the results. 
Inserting the uniform expansion in 1)45(1 we get 

+ wl^i) + ("^ + + 2)[v;°+i + Mn^„l+i + v^n+i) + A'nv;f,+i + vl^^}] + (B1) 

Note the appearance of the seemingly "higher-order" term jl^°'W^j^i. This term is, in fact, of order and therefore 
has to be included. Terms of this kind appear in other equations and the boundary conditions and one has to be 
aware (and cautious!) of their presence. The continuity equation splits into, 

lyo'+i = 0, at 0(A-") (B2) 
r/i"VF^Vi + ^^+1 = (m+l)(m + 2)1/^+1, at 0(1) (B3) 

rPF^Vi + M^^+i + T^„\+i = (m+l)(m + 2)(y„\+i + V;,Vi)-rM°W^2Vi, &t 0{fi-) (B4) 
The first of these equations immediately leads to 

W^+i = (B5) 
As (|B4ll contains both "boundary layer" and background variables, we split it further into the two equations 

rW^+^ + W^_,^ = (m + l)(m + 2)y„Vi (B6) 

rfi''W^+, = (m + l)(m + 2)V^^^ (B7) 

The same procedure will be applied in several equations to follow. Note that it should not be applied in the case of 
boundary conditions, where one must keep "boundary layer" and core variables together. At 0{jl^"') the resulting 
background continuity equation is, 

rW^+i + = {ni + l)(m + 2)V,l+, (B8) 

The respective "boundary layer" equation will not be needed in our calculation. 

Analogously we find that the radial vorticity equation implies ((28(1 . with cr — > uo, at leading order. To derive the 
next order equation (which should govern the leading order flow) we need to determine the exponent a. Potential 
terms for this equation are, 

fi{m + 1)U°;' + |ao(Ti(m + 1) + ^Mm(m + i f"^ ((7° + C/° ) 

-2aor!Q™+i(m + 2)V;^,+i + (m + l)/i'+"C/,^," = (B9) 

After some experimentation we find that a — 1/2 is the appropriate choice (other possibilities do not lead to consistent 
solutions). Then we get at 0(1); 

- ilim + 1)U^' + 2naoQrn+iim + 2)V,l^, = (BIO) 

while at ©(/i^/^) we have 

^ flim + l)C/i„ + 2f7aoQ,„+i [(m + 2)V,i+, + W^„\+i] = (m + l)fi-^/^aocJiU°, 

-217(71/2-1/2 (m + 2)Q„,+iV;0+, (Bll) 

2nii'^^Q,„+,[W^+, + {m + 2)V^,^,] - {m + l)a,U^ (B12) 
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One level up, at 0{fl), we have 

/2"i/2(m + l){al + aoa2)U^ + (m + l)aiaoU^ = -jji^'^im + 1) 
+2r!aig„+i [W^^+i + (m + 2)F^+i] + 2naoQ„^+ifi''^ [W^^+i + (m + 2)F^+i] (B13) 



and 

P^'^{m + l){/2;' = 2r!aog,„+iA'/'[W^™+i + (rn + 2)V;2^i] + 2r!aig„+i(m + 2)1/^+1 - (m + 1)^0^1 (B14) 

In the case of the "divergence" equation the leading order background flow is governed by the standard result H29|l , 
with Urn -> U^m- The 0{il~^/'^), 0{fl°) and 0{fl^^^) equations are, 

(m + l)(m + 2)[AV^'^i + a^K^'+i] - ^^<JornV^+i ~ 2f7aom(m + 2)0„,+i[/°' - (B15) 
(m + l)(m + 2)[rA3/2^u.^ ^ ao'^A'^'V;V+i + 2ao<y,rV,''„'+, + 2mf^aoV;°+i] - 2f^aomrA^/' [K^'+i + T4^„\'+i] 
-20aimrV;^Vi " 2naom{m + 2)Q^+i[rfi^/^U^ - (m + l)C/° ] - 2naim{m + 2)Qm+irUZ = (B16) 

and 

(m + l)(m + 2)al[rV,]:^, - W^^,] + 2f7aom(m + l)(m + 2)y„Vi = 2r!aomr[C+i + 

+2r!(Tom(m + 2)Q,„+i [rC/^: - (m + 1)1/^] + 2r!ai/i-i/2^(m + 2)Q„,+i[rC/^' - (m + l)C/° ] (B17) 

+2f7aimrV;Y+i + 2f^m(m + 2)Q,„+i [air[/i„' - ao(m + l)c7i„] - 2f^m(m + l)(m + 2)ctoF^+i (B18) 
respectively. Finally, the 0{ji) equation is, 

m(m + 2)Q^+,[rUZ - (m + \)Ul] = ^{m + l)(m + 2)[rV;V+i - + 2rV^+i 

+m(m + l)(m + 2)y^+i - [mrW^^Vi + + 2)M^™+i] (B19) 

To complete the required set of equations, we consider the ^-component of the Euler equation which leads to the 
pressure eigenfunction (|^ at leading order. Meanwhile the C(/i°) and ©(/i^/^) equations give 



~Kn+i = [-mao{ao-217Q2^+J(7° -m/i;7°'' + ^oQm+i{(rn + 2)ao + 2ml)}V;^+i 

Oil 



and 



+ p^{m + 2)Q„,+^V:;:'+^] (B20) 



h]n+i = T—o\n [-'^^^^lr^'^Wo-^Ql^+M.~mao{<Jo-2nQ^,^+^}Ul 

(m + Z)Lj,n+l 

+ (ToQm+i{^To(m + 2) + 2TOl]}K,\+i] (B21) 



1 

(m + 2)Qrn+l 



hln+1 = 7 — [-mflU^' + {m + 2)Qrn+iiiV^[i - rnao{ao - 2nQl^^i}U^ 



^oQm+iWaim + 2) + 2mn}Vi+^ - 2maifl-^/^{ao - f^g^+jC/^ 

2aiil-'^^Qm+i{im + 2)ao + mn}V^^^] (B22) 
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The 0{il) equations are, 
1 



"-m+l 



and 



"'m+l 



(m + 2)Qm+i 



-mao{ao - 2nQl^_^_^}U^ ^"^^^ 



{ao{m + 2) + 2mn}V,l^, + 



Qn.+iWoim + 2) + 2mn}Vl^^ 



(B23) 



1 



(m + 2)Qm+i 



2ma 



mao{ao - 2nQl^^^}Ul - ^Ty^i^o - QQl^^^jU^^ + (7oQ,n+i{Mm + 2) + 2mQ}Vl+^ 



+ ^Qrn+i{Mm + 2) + 2mr!}V;i+i - mfiU^' + fi{m + 2)Q„,+iV^" 



m+l 



(B24) 



APPENDIX C: BOUNDARY CONDITIONS 



In this Appendix we expand the boundary conditions (|22|) - l|26|) with respect to jl. 

From the conditions pertaining to the tangential traction at both r = Rc and r ~ R we obtain, 



yO/ 



l/2yV 



, at Oifl-^/^) 
2{Ul + C/O ) , at OijP) 
2(C/^ + c7iJ , at 0(^1/2) 
, at 0(/i-i/2) 



= 2v: 





m+l ' 



at Oiji'') 



at 



Continuity of the radial displacement at r = i?c demands. 

Note that the first of these conditions is trivially satisfied since we are assuming that the leading order flow 
axial. 

Finally, from the radial traction conditions at the crust-core interface and the surface we get, 



''m+l ^ ^m+lK-fi^c 



tic 



,and 



7;° + 7;° 



m+l 
"•m+l 



m+l 
"•m+l 



= -|r:+l 



"m+l ^ "m+l 



R^ ' 



W^m+l] 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 
(C6) 

(C7) 
(C8) 
(C9) 

is purely 

(CIO) 
(Cll) 

(C12) 

(C13) 
(C14) 

(C15) 



respectively. 



19 



ACKNOWLEDGMENTS 



This work was supported by PPARC through grant number PPA/G/S/2002/00038. NA also acknowledges support 
from PPARC via Senior Research Fellowship no PP/C505791/1. 



[1] N. Andersson, ApJ 502, 708 (1998) 

[2] J.L. Friedman and S.M. Morsink, ApJ 502, 714 (1998) 

[3] N. Andersson and K.D. Kokkotas, Int. J. Mod. Phys. DIO, 381 (2001) 

[4] K. Glampedakis and N. Andersson, MNRAS 371, 1311 (2006) 

[5] P. Haensel, K.P. Levenfish and D.G. Yakovlev, Astron. Astrophys. 381, 1080 (2002) 

[6] A. Rciscnoggor and A. Bonacic, Phys. Rev. Lett 91, 201103 (2003) 

[7] M. Nayyar and B.J. Owen, Phys. Rev. D 73, 084001 (2006) 

[8] N. Andersson, G.L. Comer and K. Glampedakis, Nucl. Phys. A 763, 212 (2005) 

[9] U. Lee and T.E. Strohmayer, Astron. Astrophys. 311, 155 (1996) 

[10] P.N. McDermott, H.M. Van Horn and C.J. Hansen, ApJ 325, 725 (1988) 

[11] T.E. Strohmayer, ApJ, 372, 573 (1991) 

[12] S. Yoshida and U. Lee, ApJ 546, 1121 (2001) 

[13] Y. Levin and G. IJshomirsky, MNRAS, 324, 917 (2001) 

[14] F. Douchin and P. HacnscI, Astron. Astropliys. 380, 151 (2001) 

[15] K.H. Lockitch and J.L. Friedman, ApJ 521, 764 (1999) 

[16] L.D. Landau and E.M. Lifshitz, Theory of Elasticity, Course of Theoretical Physics Vol. 7, 3rd edition (Butterworth, 
Heinemann, Oxford 1999). 

[17] CM. Bonder and S.A. Orszag, Advancx.d Mathematical Methods for Scientists and Engineers (McGraw-Hill International 

Editions, 1978) 

[18] J. Kinney and G. Mcndcll, Phys. Rev. D 67 024032 (2003) 

[19] G. Mendell, Phys. Rev. D 64, 044009 (2001) 



